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Abstract 

Injection  molding  of  semi-crystalline  polymers  has  been  one  of  the  most  important 
yet  complicated  fabrication  processes  that  combines  polymer  rheology,  heat  transfer, 
and  crystallization  kinetics.  3D  mathematical  modeling  of  polymer  flow,  and 
computer  simulation  of  distributions  of  crystallinity  in  a semi-crystalline  polymer  in 
injection  molding  applications  were  studied  in  this  project.  Our  modeling  includes 
the  shear  stress  and  thermal  conduction  in  the  width  direction  and  will  also 
eventually  provide  details  of  the  polymer  flow  in  intricate  small  cavities  of  changing 
cross  sections.  The  simulation  developed  in  this  study  predicted  3D  crystallization 
properties  of  syndiotactic  polystyrene  as  a function  of  position  in  an  injection 
molded  part.  The  predicted  crystallinity  distributions  were  effectively  described 
across  a profile  at  various  locations  along  the  flex  bar  as  a function  of  different 
molding  temperatures  and  different  hold  times. 

1.  Introduction 

Computer  simulations  of  injection  molding  of  semicrystalline  polymer  components 
have  been  an  active  area  of  research  n'6).  The  primary  goal  of  this  project  was  to 
develop  3D  computer  simulated  virtual  integrated  prototyping  (VIP),  rather  than 
traditional  2D  Hele-Shaw  approximations  or  smaller  scale  physical  trials,  to  enhance 
downstream  plastics  manufacturing-particularly  injection  molding.  Primary  semi- 
crystalline polymers  evaluated  include  polyphyenelene  sulfide  and  syndiotactic 
polystyrene.  Preliminary  results  from  this  study  have  been  presented  and  published 
elsewhere  (7  9).  Several  other  authors  have  also  published  simulation  results'10921. 
The  simulation  results  in  this  study  have  focused  almost  exclusively  on  syndiotactic 
polystyrene  (244,000  Mol.  Wt.).  The  confirmation  of  our  model  also  involved  the 
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evaluation  of  injection-molded  parts  made  from  syndiotactic  polystyrene.  The 
predicted  crystallinity  distributions  were  compared  with  measured  values  from 
molded  parts. 


2.  Governing  Equations 


Let  L,  W,  H be  the  length,  width  and  thickness  of  a slit  mold  region,  respectively. 
We  consider  the  following  mass,  momentum  and  energy  balance  equations: 
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where  the  heat  capacity,  Cp , is 

cp-c„-i*jf£ 

and  u,  P,  T,  k,  p,  Q , , X , jj,  y are  velocity,  pressure,  temperature,  thermal 
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conductivity,  density,  specific  heat,  heat  of  crystallization,  fraction  of  crystallinity, 
viscosity  and  shear  rate,  respectively. 

Many  applications  in  injection  molding  require  processing  polymers  at 
temperatures  below  the  glass  transition  temperature,  T„,  and  at  temperatures  over  the 
melting  point,  Tm.  We  propose  a new  model* l5,14>  that  allows  appropriate 
crystallization  results  to  be  obtained  over  the  entire  temperature  range  of  processing: 

Log10(-Log10/c(7) ) = Logi0(-Logi0kp)  + aY+  bZ+cV3  +01/  + e/+  fV6  (5) 

where  Y=  (T-TP)2,  k (7)  is  the  crystallization  rate  at  processing  temperature  T,  kP  is 
the  peak  crystallization  rate  constant,  a,  b,  c,  d,  e and  f are  numeric  constants  and 
Tp  is  the  peak  crystallization  temperature. 

The  sensitivity  of  kP  and  7p,  to  shear  stress,  x,  is  shown  by  the  following051: 

Tp  = Tpq  + Esik(x),  Log  io  kP  = Logio  kn+  ESk(x)  (6) 

Since  the  shear  stress  constants  ESk  and  ESTk  for  syndiotactic  polystyrene  were  not 
available  for  this  study,  they  were  assumed  to  be  comparable  to  those  by  Hsiung(13) 
for  polyphenylene  sulfide  based  on  the  data  of  Haas.et  al<16).  On  the  other  hand,  the 
true  induction  time  evaluated  must  be  generated  from  an  equation  introduced  by 
Sifleet  et  al(17). 

In  our  study,  the  induction  time  was  determined  using  the  following  model: 
ti  = tib  e0,i'T  ~Tb^ 


(7) 
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The  effect  of  shear  stress,  x , on  the  induction  time  was  again  addressed  from  the 
expressions  proposed  by  Hsuing'151  for  polyphenylene  sulfide  as: 

Tb  = Tbq  + EsTii(T),  Logiotib  = logio  tibq  - Esti(t)  (8) 

Values  of  the  constants  for  ESTti  and  ESti  were  evaluated  from  the  data  of  Hsuing051 
and  related  studies06’18  20)  for  polyphenylene  sulfide. 

The  melt  viscosity,  t|,  of  the  syndiotactic  polystyrene  was  evaluated  using 
the  Cross  law<2l)  as  a function  of  shear  rate  (see  Verhoyen  et  al (22)). 


2.1.  Fractional  Crystallinity  System  of  Integro-Differential  Equations 


It  was  known(8'9)  that  the  following  system  was  satisfied  by  the  fractional 
crystallinity, 
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I = J K (T  (x,  y,t),  T(x,  y,t))dt 

rcx 

Initial  polymer  melt  temperature  T0  was  taken  as  the  polymer  processing  temperature 
at  the  shutoff  nozzle.  It  was  observed  that  the  degradation  of  syndiotactic  polystyrene 
is  considerably  lower  if  the  temperatures  were  kept  below  305°C.  The  melt 
temperature  was  taken  to  be  305°C.  T0  = 305°C  while  wall  temperature,  Tw,  is  a 
constant. 


T(x,  H,  z,  t)  = T(x,  0,  z,  t)  = T (x,  y,  W,  t)  = T(x,  y,  0,  t)  = 7W  (11) 

The  following  homogeneous  velocity  on  the  boundaries  of  the  mold  was  considered. 

u(x,H,  z,t)  = u(x,0,z,t)  = u(x,y,  W,t)  = u(x,y,0,t)  = 0.0  (12) 

The  melt  front  of  the  polymer  was  considered  to  be  flat  and  the  temperature  was 
assumed  to  be  uniform  all  over  the  melt  front.  The  uniform  temperature  is  equal  to 
the  centerline  temperature  at  the  stream-wise  location  immediately  upstream  of  the 
melt  front.  The  curvature  and  the  transverse  flows  associated  with  the  melt  front  due 
to  the  fountain  flow  effect  were  neglected.  Thus  we  have  the  following  equality: 

T(xmf>y>z>U-T^xm,  - Ax,  ^ yy  > 1 j U3) 


2.2.  Numerical  Method 


640 


R.  D.  SUDDUTH,  P.  K.YARALA,  Q.  SHENG,  K.  NICHOLS 


Introduction  of  the  discrete  parameters.  A*,  Ay,  A z and  At  along  x,  y,  z (i,  j,  k) 
direction  s and  time,  t(n),  into  the  energy  balance  equation  yields  the  following 
implicit  schemes  via  the  parallel  splitting  algorithm'8,9,235. 
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Thus,  by  applying  the  parallel  splitting  algorithm'235,  the  solution  of  the  energy 
balance  equation  can  be  obtained  with  a neglectful  error  as 


T n+1  = 1 / 2(exp{  AtM ! } exp{  AtM  2 } exp{  AtM  3 } + 


exp{  AtM3}  exp{  AtM2}  exp{  AtM^T  n (14) 

where  M|,  M2  and  M3  are  tri-diagonal  matrices  generated. 


4.  Computer  Program  Simulation  Results 

4.1.  Simulation  of  Crystallinity  Distributions  at  Different  Holding  Times  and 
Different  Locations 

Computer  simulated  crystallinity  distributions  at  a mold  temperature  of  200°C  and 
an  injection  speed  of  23.2  cc/sec  at  Location  #1  in  Figure  1 are  summarized  in 
Figures  2 as  a function  of  holding  time.  In  all  the  figures  plotted  from  the  simulation 
results  the  direction  shown  as  Z direction  corresponds  to  width  (1/2  in)  of  the  sample 
and  the  direction  shown  as  Y direction  corresponds  to  thickness  (1/8  in)  of  the 
sample.  The  location  indicated  as  ‘entrance’  is  located  near  the  gate  and  the  location 
indicated,  as  Location  #3  is  located  further  down  the  sample  along  the  direction  of 
flow.  The  third  axis  in  the  figures  shows  the  magnitude  of  the  fraction  of 
crystallinity,  X,  developed.  The  value  of  X = 1 along  this  axis  direction  indicates 


MEASUREMENT  AND  3D  COMPUTER  SIMULATION. . . 


641 


the  maximum  crystallinity  that  can  develop  at  any  point.  The  maximum 
corresponded  to  a heat  of  fusion  which  was  approximately  61.7  percent  of  the 
maximum  crystallinity  possible  or  fc  = .637  giving:  |fc  A^|  = .617  (53.2 
Joules/gram)  = 32.82  Joules/gram. 

As  shown  elsewhere113’14',  the  maximum  crystallization  rate  for  the 
syndiotactic  polystyrene  utilized  in  this  study  can  be  achieved  at  a temperature  of 
approximately  200°C.  A rapid  surface  crystallinity  development  at  200°C  is 
indicated  in  Figure  2 as  a function  of  holding  time  at  location  #1 . The  computer 
simulated  crystallinity  at  a mold  temperature  of  200°C  and  an  injection  speed  of  23.2 
cc/sec  is  also  indicated  at  several  different  locations  in  Figure  3 at  a holding  time  of 
60  sec. 


4.2.  Crystallinity  Distribution  at  Two  Locations  as  a Function  of  Mold 
Temperature  at  a Given  Holding  Time  and  a Specific  Injection  Speed 


The  injection  molding  we  considered  was  at  mold  temperatures  of  50°C,  90°C, 
150°C  and  190°C  and  an  injection  speed  of  23.2  cc/sec.  The  simulations  for  these 
runs  are  indicated  in  Figure  4 for  location  position  indicated  in  Figure  1 as  the 
entrance.  Assuming  that  the  mass  forming  the  injection  molded  polymer  bar 
remained  the  same  independent  of  the  temperature  of  the  mold,  then  the  old  mold 
temperature  and  new  estimated  mold  temperatures  were  as  follows: 

Tmold.°C  TNEw,°C 

50°C  132.3  (135) 

90°C  159.4(160) 

150°C  200(200) 

190°C  227.1(230) 


The  new  computer  injection  molding  simulations  for  mold  temperatures  of  135°C, 
160°C,  200°C  and  230°C  are  indicated  in  Figures  5 at  the  entrance  location. 


5.  Experimental 

The  syndiotactic  polystyrene  used  in  this  study  was  supplied  by  the  Dow  Chemical 
Company  and  had  a molecular  weight  currently  characteristic  of  the  production 
grade  of  syndiotactic  polystyrene  (Mol.Wt  = 244,000).  Most  of  the  experimental 
details  for  the  crystallinity  and  injection  molding  measurements  have  been  reported 
elsewhere03'14’.  The  crystallinity  measurements  determined  at  four  different  mold 
temperatures  (50°C,  90°C,  150°C  and  190°C)  are  summarized  in  Figures  6 and  7. 

Since  the  maximum  crystallization  rate  was  found  to  occur  at  approximately 
200°C,  it  was  expected  that  the  mold  temperature  of  190°C  would  have  given  the 
maximum  amount  of  surface  crystallinity.  However,  the  maximum  surface 
crystallinity  was  achieved  with  a mold  temperature  of  150°C  instead  of  190°C.  This 
result  strongly  suggested  that  the  effective  mold  temperature  at  150°C  was  probably 
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nearer  to  the  200°C  as  proposed  in  an  earlier  section.  Comparison  of  these 
experimental  results  with  the  simulated  results  in  Figures  4 & 5 would  also  be 
consistent  with  this  observation. 

6.  Conclusion 

In  this  study,  four  mold  temperatures  (50°C,  90°C,  150°C  and  190°C)  were 
experimentally  evaluated  and  compared  with  computer-simulated  results.  For  the 
injection  molded  bar  simulated  in  this  study  it  was  found  that  the  maximum 
crystallinity  characteristic  of  a theoretical  mold  temperature  of  approximately  200°C 
was  experimentally  achieved  at  a mold  temperature  of  approximately  150°C  at  an 
injection  speed  of  23.3  cc/min.  Using  a correction  factor  approach  and  local  heating 
it  was  found  that  the  thermal  transfer  of  energy  apparently  allows  a temperature  of 
approximately  200°C  to  be  achieved  in  the  mold  at  a mold  temperature  of  150°C. 

By  addressing  effective  computer  simulation  and  appropriate  experimental 
outputs  for  syndiotactic  polystyrene,  it  has  been  shown  that  an  optimized  control  of 
the  injection  molding  process  can  be  developed  for  better  commercial  applications. 
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Figure  3 Crystallinity  Distribution  at  Different  Locations  at  a 200°C  Mold 
Temperature  after  a Hold  Time  of  60  Seconds 


Figure  4 Crystallinity  Distributions  at  the  Entrance  Location  at  Four 
Different  Mold  Temperatures  (50°C,  90°C,  150°C,  190°C)(Injection  Speed  of  23.2 
cc/sec  and  a Hold  Time  of  60  sec) 
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Figure  5 Crystallinity  Distributions  at  the  Entrance  Location  at  Four 
Effective  Mold  Temperatures  (135°C,  160°C,  200°C,  230°C)  (Injection  Speed  of  23.2 
cc/sec  and  a Hold  Time  of  60  sec) 


Fraction  of  Dlslanca  from  Surfaco  Relative  to  ThicfoMS 


Figure  6 Percent  Crystallinity  vs  Distance  from  the  Surface  at  Entrance  Location  for  Three  Mold 
Temperatures  (90aC,  150aC,  190aC) 


Fraclon  of  Ostanca  From  Surfaco  Ratafivo  to  Thickness 


Figure  7 Percent  Crystallinity  vs  Distance  from  the  Surface  at  Entrance  Location  and  Two  Mold 
Temperatures  (50°C,  90°C) 


